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Abstract. We investigate a large class of antiferromagnetic spin-1 chains with nearest neighbour interaction 
and exactly known matrix product ground state. The spectrum of low-lying excitations is calculated 
numerically by DMRG and exact diagonalisation. Spin liquid behaviour with an excitation gap is observed 
meeting the Haldane scenario. Further, we compare properties of the anisotropic model with the well- 
known isotropic AKLT model and use the analytical single mode approximation to obtain a quantitative 
understanding of the excitation gap. The low-lying excited states can be interpreted in terms of a magnon- 
like elementary excitation. 

PACS. 75.10.Jm Quantized spin models. 



1 Introduction 

Quantum spin systems have been studied intensively in 
recent years (for a review, see e.g. ^). Special interest 
has been taken in the case of low dimensions where strong 
quantum fluctuations destroy magnetic order j2j . Antifer- 
romagnetic spin chains with S = 1 have been the focus of 
theoretical and experimental studies since Haldane's con- 
jecture |3| about the difference between integer and half- 
integer spins. Antiferromagnetic integer spin chains are 
predicted to be generically massive while half-integer spin 
chains have a vanishing gap which can be rigorously shown 
for the isotropic spin-^ chain using the Bethe ansatz so- 
lution [3]. For the isotropic Heisenberg S" = 1 chain with 
nearest neighbour exchange no exact solution has been 
found, but there is convincing evidence from numerical 
studies that a finite gap exists in the thermodynamic limit 



In a most general spin-1 Hamiltonian, bilinear as well 
as biquadratic interaction terms along with possible an- 
isotropics emerge. Consideration of only bilinear exchange 
with anisotropic interaction and single-ion anisotropics 
leads to a well-known two-parameter phase diagram which 
includes massive as well as massless phases (see e.g. [H] and 
references therein). The spin-1 Hamiltonian with isotropic 
nearest neighbour - bilinear and biquadratic - interaction 
has also been focus to extensive theoretical studies, see 

including the Heisenberg Hamiltonian exhibits a variety 
of magnetic phenomena like ferromagnetism, dimerization 
and a gapped spin liquid (Haldane) phase. A valence- 
bond-solid model (AKLT model) with exactly known ground 
state representing this phase has been introduced by Af- 



fleck et al. 'TO. Only a few isotropic models with antifer- 
romagnetic bilinear and biquadratic exchange possess ex- 
actly known spectra |17II18II19II21)II21| . While experiments 
on many spin-1 compounds |22ll2;-{ll24ll25ll2f)j suggest that 
biquadratic exchange can usually assumed to be small, 
Mila and Zhang (221 presented a microscopic model with 
non-negligible biquadratic exchange in order to under- 
stand a massless phase of the spin-1 vanadium oxide LiVGe206- 

In I2B] a large class of spin-1 models possessing a unique 
matrix product ground state was introduced. This model 
class includes biquadratic terms and anisotropic exchange 
at the same time and has been proven rigorously to ex- 
hibit exponentially decaying spin-spin-correlations. In the 
isotropic limit it reduces to the AKLT model jl^. While 
the ground state and some of its properties have been cal- 
culated exactly the excitation spectrum in the anisotropic 
case has not been subject of studies yet. It has only been 
known that generically a finite gap is observed, leading to 
spin-liquid behaviour. 

In this work we present numerical and analytical evi- 
dence for a finite excitation gap and discuss the spectrum 
and structure of low-lying excitations. Results from an an- 
alytical single mode approximation capture the nature of 
the lowest excited states and give an exact upper bound 
to the excitation gap. 



2 Model and symmetries 

We consider a spin-1 chain with rotational invariance in 
the {x, y)-plane, invariance under —^ —S^ and trans- 
lation and parity invariance, i.e. invariance under the ex- 
change J <-> J + 1. The most general Hamiltonian with 
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nearest neighbour interactions satisfying the above sym- 
metries 



+ J4 [{S^S^+, + SjSjVi) + h.c' 

+ D {Sf + S]^,') + c} 



(1) 



gives rise to a model class with seven parameters including 
an energy off-set c and a scale, reducing the number of rel- 
evant parameters to five. The additive constant c is used 
to shift the ground state energy to Eq = 0. The Hamil- 
tonian |^ includes the bilinear, anisotropic spin-1 chain 
(J2 = J3 = J4 = 0) as well as a chain with only isotro- 
pic bilinear and biquadratic exchange {A = 1, J2 = J3 = 
J4,D — 0), both comprising the Heisenberg model. 

A large subclass of the general model class has been 
found [2H| to possess a unique ground state which can be 
written as a product 



|Vo) =tr(77i(i) 



(2) 



of state matrices 



m 



U) 



a 7^ 0, cr ± 1 (3) 



characterised by the real parameter ^ a and the discrete 
parameter ct = ±1. Minimizing not only the total energy, 
but the energy of the local two-site Hamiltonian as well, 
the matrix product ground state (MPG) belongs to the 
class of optimum ground states [2H|- The parameter a 
represents a degree of freedom connected to the aniso- 
tropy of the model. Further, taking excited states into 
account, three positive real spectral parameters Aq, Ai, A2 
complete the parameter space of the MPG model. Neglect- 
ing a trival scale three relevant real parameters and one 
discrete parameter survive. For models with MPG a con- 
venient representation of the spin-1 Hamiltonian makes 
use of spectral parameters and projection operators (see 
appendix) rather than spin operators Q and shall be em- 
ployed in this paper. 

Referring to the spin-1 model class, among others, two 
unitary transformations 



p2 



MS' 



(4) 



with A denoting the subchain of even sites, mapping the 
parameter space onto itself can be found. Corresponding 
Hamiltonians show identical energy spectra, thus reduc- 
ing the parameter space we need to discuss. The unitary 
transformation 7?^ maps the parameter a on the model 



^ Note that we choose the parameter o in a different way 
compared to |28| in order to stress the spin-flip symmetry of 
the state matrix m^-'\ 



with —a. It has no impact on the momentum. Therefore, 
in the following we concentrate on \a\. The second uni- 
tary transformation changes model parameters (a, a) 
to {—a,— a). Please note that the momentum of odd-S'^ 
states is modified from k to tt — k whereas the momen- 
tum of states with even remains unchanged under this 
transformation. 

While the ground state is determined by only the an- 
isotropy parameter a and the discrete parameter cr, exami- 
nation of the excitation spectrum requires consideration of 
the spectral parameters as well. Choosing the parameters^ 
Ao = Al = A2 = 1, a = —2, a = —1 , an isotropic point of 
the MPG model class is reached. This special case corre- 
sponds to the AKLT model and has been studied inten- 
sively Therefore, with fixed Aq = Ai = A2 = 1 
we vary the parameter \a\. 



3 Numerical approach to low-lying excitations 

Low-lying excitations of the model are obtained by us- 
ing density matrix renormalisation group (DMRG) |31II32| 
and exact diagonalisation data. Our DMRG-program im- 
plements the infinite system algorithm with conservation 
of the quantum number as well as periodic boundary 
conditions. Though usually a higher accuracy of DMRG 
data is achieved for open boundary conditions we have 
chosen periodic boundary conditions for a better compar- 
ison with analytical calculations for spin chains of different 
lengths. Moreover, dealing with a gapped system includ- 
ing a unique matrix product ground state the calculated 
energies are supposed to converge quickly This has 
already been observed for the isotropic case in the vicinity 
of the AKLT point |lllll2j . In order to check the accuracy 
of our data we tested the convergence of the numerical 
method related to the chain length L and to the number 
TO of states kept in each DMRG-iteration as well. With 
L = 100, m < 120 reliable gap information for the MPG 
model class could be achieved. Since in real space DMRG 
momentum is not preserved we used exact diagonalisation 
of smaller chains {L < 14) to discuss the momentum of 
the excitation levels. 



3.1 Isotropic interactions (AKLT model) 

Before we discuss the general case we briefiy recall the 
results for the isotropic Hamiltonian (a = —2, a — — 1). 
In fact the excitation spectrum of the AKLT model has 
been studied using numerical and analytical methods |12l 
I34| and is well understood. The elementary excitations 
are magnon-type spin-1 particles. The excitation gap is 
determined by the lowest energy Es{k) of the triplet exci- 
tations {S — I) with momentum k = n which belongs to 
a discrete branch below a continuum. The second excited 
state at fc = TT has energy 3i5i(7r) and spin 5 = 3. On the 
other hand, the gap at momentum A: = is twice as large 
as that at A: = TT and has spin S = 2 ^21- This suggests 

^ This equals Ji = |, Zi = 1, D = 0, J2 = J3 = J4 = 
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Fig. 1. Lowest excitations in the = (o), 1 (a), 2 (□) sub- 
spaces plotted versus anisotropy parameter \a\ (DMRG-data). 



that the corresponding states are scattering states of three 
resp. two magnons with S = 1 and momentum k = tt. 



3.2 Influence of anisotropy 

In fig. n the numerically calculated lowest excitations in 
three 5'^-subspaces are plotted as function of the param- 
eter \a\. For |a| 7^ a finite excitation gap to the unique 
ground state is observed |2HI ■ At the point a — the model 
reaches a higher symmetry (invariance under transforma- 
tion V^) leading to a three- fold degenerate ground state 
with vanishing excitation gap in the thermodynamic limit. 
At the isotropic AKLT model with a = —2 and a = —1 
the lowest excitation is a triplet with total spin 5=1 and 
momentum fc = tt, in agreement with previous results (see 
Sec. l3.1|l . A degeneracy of the energy levels in the = 0,1 
subspaces can thus be seen at \a\ = 2. For smaller val- 
ues of |a| the anisotropy is of easy-plane type leading to 
a splitting of the lowest excitation level. Here, excitations 
along the z-axis are energetically favourable. Therefore the 
lowest excitations are found in the — 1 subspaces. In 
contrast, for \a\ > 2 an easy- axis anisotropy favours exci- 
tations with a vanishing total S^. Further, we observe a 
level crossing in the — subspace at a = a; « 1.4. An 
obvious conjecture, consistent with the numerics, would 
be ai = \f2. At this point the interaction parameter J3 
changes sign. However, for a different choice of anisotropy 
parameters Aj the level crossing seems not to be related 
to such a change of sign. It can even be observed in cases 
where there is no sign change at all (see Sec. I3.4f) . 

For smaller \a\ the lowest excitation with vanishing S^' 
degenerates with the lowest energy levels in the S*^ = 2 
subspaces. Its energy is close to double the size of the = 
1 excitation gap. This gives rise to the conjecture that the 
lowest excitations in the 5^ = 0, 2 subspaces consist of two 

= \ excitations. In the following sections we find more 
evidence for this picture. For a > a/ a different energy level 
in the 5"^ = subspace becomes energetically favourable 
which, for a — s- cx), decreases to zero. 



With DMRG, we calculated up to five lowest eigenen- 
ergies of the system in each S'^-subspace. All higher exci- 
tation energies were found to lie, in the thermodynamical 
limit, infinitesimally close to the lowest excitation energy 
in each subspace in accordance with an analytical argu- 
ment of Fath and Solyom Jl] . 



3.3 Dispersion 

As a method in real space, DMRG gives no information 
about the momentum of the excitations. In order to have a 
closer look at the dispersion curve, the low-energy spectra 
of smaller chains {L < 14) were calculated by exact diag- 
onalisation. Figures |21 and |21 show the dispersion of low- 
lying excitations at four representative points for a = —1 
(Fig. [SJ and a = 1 (Fig. O. In order to estimate the ac- 
curacy of the exact diagonalization results we have also 
included the DMRG data for the energy gaps. They agree 
very well with extrapolations from results for systems with 
L < 14. Interestingly, large finite-size corrections are only 
found in the 5^ = subspace for small |o|. Later we will 
argue that the corresponding states can be interpreted as 
scattering states of two magnons. The other states found 
using DMRG are single magnon excitations. Their ener- 
gies converge evidently faster with increasing system size. 

In the following subsections we will discuss the spectra 
for the regions < \a\ < ai, ai < \a\ < 2 and \a\ > 2 
separately. Figs. [21 and El show typical examples for each 
region. 



3.3.1 < |a| < ai 

For anisotropy parameters < |a| < a; the lowest exci- 
tation is found in the subspace of S"^ = 1 and has mo- 
mentum k — TT {a — —1) or k ^ {a ~ 1). The lowest 
states in the S"^ = and = 2 subspaces are degener- 
ate. They have momentum k — and an energy which is 
twice the value of the gap. The corresponding states might 
therefore by interpretated as scattering states of two ele- 
mentary magnon- type excitations with = 1. 



3.3.2 ai < |a| < 2 

As already mentioned in Sec. 18.21 at a; « 1.4 a level cross- 
ing is observed for the second excited state. In the re- 
gion a; < \a\ < 2 the scattering states for small momenta 
still exist. However, another state with = has now a 
lower energy. The momentum of the corresponding second 
gap, which is smaller than twice the gap in the = 1- 
subspace, for cr = —1 is now given by = tt instead of 
k = 0, and vice versa for cr = 1. Later (in Sec. 14. 1|) we 
will show that this new state is well described by a single 
mode approximation and therefore corresponds to another 
single-magnon excitation. This state is result of the split- 
ting of the Haldane triplet due to the anisotropy. 
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Fig. 2. Dispersion of lowest excitations in the subspaces for = (o), 5^ = 1 a and 5^ = 2 (□) at a = 0.5 (top left), 
o = 1.5 (top right), a = 2 (bottom left), a = 5 (bottom right) for u = —1. The lowest excitation energies in each subspace are 
indicated by a dotted line which only serves as a guide to the eye. The data are obtained by exact diagonalisation {L = 14). 
For comparison the DMRG results (x) from Fig. have also been included. 



3.3.3 \a\ = 2 

For the \a\ = 2, cr = —1, including the isotropic point for 
a = — 2, the well-known properties of the AKLT model 
are reproduced (Sec. l3.l|l . The model's isotropy leads to a 
degeneracy of the = 0, 1 excitation branches indicating 
the triplet character of an excitation with total spin 1. As 
can easily be seen, the gap still occurs at a momentum 

k — TT. 

For |a| = 2 and cr = 1 the model is not isotropic and 
hence no degeneracy of the = 0, 1 states is observed. 
However, the transformation which changes a — > —a 
also transforms the momenta k of odd-S*^ states to tt = A;. 
This leads to the observed symmetry of the lowest excited 
states around k = ^. 

3.3.4 \a\ > 2 

For \a\ > 2 the structure of the low- lying excitations 
changes compared to the case \a\ < 2. Now the lowest exci- 
tation is found in the S*^ = subspace. We will later show 
that it again can be interpreted as a magnon. However, 



the energy gap still occurs for momentum k = tt. Also the 
second excitation, which now is found in the — 1 sub- 
space, takes its minimum energy at k — t: (for a = — 1) 
resp. k = (for cr = 1) . Note that for a large anisotropy 
(|a| — 5 in Figs.EJ all low- lying excitations (except for 
the first excitated state) have almost no dispersion. 

3.4 Variation of the spectral parameters 

In fig. ^ DMRG data for models with varying spectral 
parameters Aq, Ai, A2 are plotted. While increasing one of 
the parameters to A^ = 2, the others are kept constant 
{Xj^i = 1). All diagrams share the qualitative appearance 
of the original model with Aq = Ai = A2 = 1 (fig. 0) 
containing the isotropic point. Again, we observe a level 
crossing in the — subspace. A second level crossing 
occurs between levels in the subspace with — and 
5'^ = 1 at the point, which represents the isotropic system 
in the original model. Note that for the increase of any 
A-parameter all excitation levels are lifted compared to 
the original model. Yet, the different excitation branches 
show different scaling behaviour under the variation of the 
spectral parameters. 
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Fig. 3. Same as Fig. 13 with a = 1. 
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In the 5^ = subspace with |a| beyond the level cross- 
ing we find a significant Ai-dependence of the decreasing 
excitation branch. The influence of the other spectral pa- 
rameters Ao and A2 is less important. The lowest exci- 
tations in the S^' ~ 1 and S^' — 2 subspace - for small 
\a\ degenerate with the S'^ ~ Q increasing branch - react 
to variations of the spectral parameters in a similar way. 
While Al has strong influence in the region with small \a\ 
and while increasing A2 has greater impact in the large 
\a\ limit, there is only a weak Ao-dependence. In models 
with small \a\ states with local spin configuration with 
no z-axis alignment dominate. Here, we have only a small 
possibility to find and (— , — ) adjacent spins which 

contribute to A2-dominated excitations. This is different 
in the large \a\ regime where spins are strongly aligned 
along the z-axis leading to the observed A2-dominance. 



4 Single mode approximation 



The similar scaling behaviour of the increasing branch- 
es in the 5^ = 0, 1, 2 subspaces suggests a common nature 
of these excitation levels. Particularly, it supports the pic- 
ture that the increasing energy branches in the 5^ = 0, 2 
subspaces are formed by two weakly interacting magnon- 
type excitations {S^ = 1). 



In order to substantiate the interpretation of the elemen- 
tary excitations as magnons, we use the single mode ap- 
proximation (SMA). The SMA is a variational approach 
to find approximate low-lying excitations by making use 
of ground state correlations. Therefore, the ground state is 
locally perturbed similar to a spin-wave. This single mode 
state has a certain momentum and, since it is orthogonal 
to the ground state, its energy is an exact upper bound 
to the excitation gap. Introducing the SMA, Arovas et al. 
|29| gave an exact upper bound to the excitation gap of 
the AKLT model. In this paper, we extend their approach 
to the general anisotropic MPG model class in order to 
illustrate the spectrum of low-lying excitations. 

With the locally acting perturbation operator 0„ we 
construct a translationally invariant state 



IV'fc) 



VI 



(5) 
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where ki'o) is the exactly known matrix product ground 
state The energy 



{ipk\H\ipk) 

E„x e*'=("-"'K'Ao|ol,0„|V'o) 



(6) 



of this variational state can be regarded as a quotient of 
expectation values of the exactly known matrix product 
ground state which can be calculated via the transfer ma- 
trix method 35 . It is an exact upper bound to the excita- 
tion gap if the variational state is orthogonal to the ground 
state. Complying with conservation of the z-component of 
total spin, in the following we investigate the subspaces 
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cal DMRG data (circles) and energy expectation values of the 
single mode approximation (dashed line A; = 0, full line k = n) 
are plotted vs. anisotropy parameter \a\. For comparison the 
thin dotted line gives the = 1 excitation branch (DMRG 
data). 



with = 1 and = separately. The behaviour of the 
gap near the boundaries of the MPG region is of special 
interest where it is expected to vanish. 



4.1 Subspace = and behaviour for large \a\ 

Considering a local application of the S'^-operator, the 
thus built state 

|Vn = 4f Ee^'^^nko) (7) 
n=l 

is orthogonal to the ground state \tpo) and has a van- 
ishing total ^^-value. In order to find the lowest excita- 
tion with — the most general SMA perturbation 
operator would include contributions of (S^)'^ — a/(l -I- a) 
which is another locally orthogonal perturbation opera- 
tor. However, a superposition of both possible operator in 
this subspace is not mandatory since one can show that 
variational states with only 5"^ as a perturbation operator 
already provide the lowest SMA energy levels. 
We find the variational energy 



2 lal 



rAo + A] 



1 



(l + |a|)2 il + \a\r 



■ cos k 



(8) 



as an upper bound to the lowest excitation in the 5*^ = 
subspace. 

In the case of |a| > 2 the excitation gap is found in 
the S'^ = subspace with a momentum n. The fc = tt 
variational energy 



2\a\ 



■An + A] 



(9) 
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tends to zero for \a\ oo. At the AKLT point (|a| = 2) 
the SMA variational energy i?AKLT — 10/27 of Arovas 
et al. 1211 is reproduced. For smaU \a\ the lowest — 
excited state has a momentum fc = 0. 

Figure [51 shows the SMA energies in comparison to 
the excitation gap calculated by DMRG. For large val- 
ues of \a\ we find the best agreement between analytical 
and numerical data in considering the variational state 
with momentum k — n. While at the point of the level 
crossing (|a| w 1.4) the SMA-energy is clearly higher than 
the DMRG data, both analytically and numerically calcu- 
lated energies become asymptotically equal in the |a| — > oo 
limit. For small \a\ the lowest excited state of this subspace 
has vanishing momentum. Here, the variational state with 
fc = is not capable to replicate the numerically calculated 
excitation branch as suitably as in the large- |a| case. This 
indicates a character of the small- |a| excitation branch 
significantly different from the one-magnon picture of the 
excitation gap. 



4.2 Subspace ^ 1 and behaviour for small \a\ 

In order to describe the lowest excitation in the ^ I 
subspace we use the most general local perturbation op- 
erator 

Ol+^ =S+[b+{2b-l)S^] (10) 

comprising a variational parameter 6 S JR. For 6 = 0, 
we have 0'^+^ oc S^S^ whereas for 6 = ^ it reduces to 
O*-^^ cx . The corresponding variational state 



(11) 



is orthogonal to the ground state and has a total S^' — 1 
and momentum fc. Since in the MPG model with a = I 
we expect the excitation gap at vanishing momentum, in 
the following, we set fc = and assume a > without loss 
of generality. Square of the norm 



rgy expe( 

{i;+\H\i;+) = 



1 



1 + a 

and energy expectation value 
1 fl 



l + b^ + ^a{l-bf 



(12) 



(i + «r 

2Aoa 



-Aia^(l-6)" + 2A2a6' 



i{l - 6) - by 



(13) 



of the variational state can be calculated via the transfer 
matrix method. In the following we choose the realistic 
spectral parameter configuration Ao, Ai, A2 = 1, including 
the AKLT model. Thus, we obtain the energy 



W{1 - by + 2ab^ +2T^[a-b{l + a)] 



(14) 



(1 + a) 1 + 62 + ia(l - 6)' 




Fig. 6. Lowest excitations in the — 1 subspace. Numerical 
DMRG data (a) and energy expectation values of the single 
mode approximation (straight line) are plotted vs. anisotropy 
parameter \a\. Thin dotted lines give the — excitation 
branches (DMRG data). The inset shows the minimizing vari- 
ational parameter &min in the same range of \a\. 



of the variational state as a function of a and 6. Minimizing 
the energy of the variational state yields the condition 



4a2 



7a + 6 



a2(a-2) 

V4a5 + 36a* + 80a^ + 97a^ + 84a + 36 
a2(a-2) 



(15) 



(inset of fig. |SJ|. 

In fig.Elthe energy of the variational state is depicted. 
A good agreement with the numerical data, particularly 
for small a, is observed. A series expansion for small a 



^=^«-^«' + 0(a3), 
E+ = la'-^-a' + \a* + 0{a% 



(16) 



predicts a decrease of the excitation gap of order a2. Here, 
in the a — > limit, the optimal choice of the local vari- 
ational operator is S'+S'^. A fit of the numerical DMRG 
data in this region 



DMRG 



0.500a2 - 0.497a'^ - O.lSga"* 



(17) 



can be well reproduced by the variational calculation up 
to the third order in a. At the isotropic point (|a| ~ 2) 
we find 6niin = 5 (S""*" as best variational operator) and 
-E-AKLT = 10/27, both in accordance with the triplet exci- 
tation picture. For large a we find 



(18) 



6„iin = 1 - 2W - + 0(a- 
" a 



E+ = 1- J- + 0{a- 
" a 
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basis states 



The energy gap and lowest excitations of the anisotropic 
spin-1 chain with a matrix product ground state have been 
subject to numerical and analytical calculations. DMRG 
studies showed an excitation gap for finite values of an- 
isotropy parameter \a\. In the isotropic case the lowest 
excitation is a triplet of momentum tt. Under anisotropic 
interactions this excitation triplet splits and gives rise to 
two separate excitation branches with =0 and = 1 
establishing the gap for large \a\ and small \a\, respec- 
tively. 

In the large \a\ regime the lowest excitation can be 
identified as a magnon of vanishing and momentum 
k = TT. For \a\ < 2 the excitation gap is found in the 
5^ = 1 subspace with momentum fc = tt (fc = 0) in case 
of a parameter choice a = —I (a = I). Both numerically 
found excitation energies can be explained quantitatively 
well by the single mode approximation. From DMRG we 
find that the = 2 and the increasing = excita- 
tion branches have twice the magnitude of the = 1 
excitation gap. Together with results from numerical in- 
vestigation of the momentum and of the scaling behaviour 
under variation of the spectral parameters this is strong 
evidence that these excitations consists of two weakly in- 
teracting S^- = 1 magnons. 

At a = a phase transition to a three-fold degener- 
ate ground state with excitation continuum occurs. In its 
neighbourhood we observed numerically and analytically 
a decrease of the excitation gap with a^. 



^±2) =|±±), 

=^(|±0)+a|0±)), 

^o") =^(l + -)-|-+)), 
-^^i) = Vlfe(°[l +-) + !-+)] +2|00», 

'^02> = (I +-> + !-+> -«|00» 

(19) 

with cr = ±1 distinguishing symmetric and antisymmet- 
ric states. Considering the projection operators V!^ := 
I "^n I orAo these basis states and paying attention to 
the spin-flip symmetry we find a full representation of the 
local Hamiltonian 



W = ^2 {V2+V-2) 

+ A+ {vt+v\) + \-^ (^r+^ii) 



(20) 



which depends on six spectral parameters A2 , , A^^ , Aqj^ , 
Xq2, Aq and the superposition parameter a. The spin-chain 
has a MPG if A2, Aj^"", A(J2 > and Af , Aj^, Aq = leading 
to a local MPG Hamiltonian 

/i,- ,+1 = A2 {T2 + V-2) + Al (Pf^ + VZ^) + AoPo+2 (21) 

with Al = Xi and Aq = Ao2- Spin and spectral parameters 
of the MPG model class are mapped onto each other by 



This work has been performed within the research program of 
the Sonderforschungsbereich 608 supported by DFG. 



6 Appendix A 



Apart from the spin operator description, the Hamiltonian 
of the matrix product ground state model can also be rep- 
resented by a sum of projection operators on local ground 
states. This spectral representation reflects the model's 
symmetries and helps to understand the free parameters 
of the MPG model. 

Due to the rotational invariance in the {x, ?/)-plane 
we can divide the 9-dimensional Hilbcrt space of a two- 
site spin-1 system into five subspaccs of conserved = 
Sf + 82- The eigenstates of the 1-dimensional subspaces 
with = 2 are fixed. Further, parity invariance leads 
to discrimination of symmetric and antisymmetric states. 
In the 2-dimensional subspaces characterised by S*^ = 1 
each subspace has one symmetric and one antisymmetric 
eigenstatc. In case of 5^ = 0, however, after identifying 
one antisymmetric state, a 2-dimcnsional symmetric sub- 
space leaves one degree of freedom which is described by 
the superposition parameter a. Finally, we find nine local 



Ji 

J4 
D 



jo-Ai, 



2 






a 


2 




1 


-a2 



2 + a^ 



Ao — Al -I- -A2, 
Ao - 2'^^!^ 
Ao + ^Ai, c = 



"'^aT' ^^=2T^^°' 
1 



(22) 



a^ + 2 



Ao- 
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